Tangled Nature: a model of evolutionary ecology. 
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We discuss a simple model of co-evolution. In order to emphasise the effect of interaction between in- 
dividuals the entire population is subjected to the same physical environment. Species are emergent 
structures and extinction, origination and diversity are entirely a consequence of co-evolutionary 
interaction between individuals. For comparison we consider both asexual and sexually reproduc- 
ing populations. We also study competition between asexual and sexual reproduction in a mixed 
population. In either case the system evolves through periods of hectic reorganisation separated by 
periods of coherent stable coexistence. 

*Author to whom correspondence should be addressed. E-mail: h.jensen@ic.ac.uk 



I. INTRODUCTION 

It is difficult in experiments and observations to bridge 
the gap bet"ween ecological time and evolutionary time 
(Pimm, 1991). Nevertheless, since Dar-win's publication 
of The Origin of Species (Dar"win, 1859) it has been gen- 
erally agreed that the intricate and complex ecologies sur- 
rounding us are the product of Natural Selection operat- 
ing on vast numbers of successive generations. We kno"w 
that the slo"w gradual effect of mutations and Natural 
Selection is the long term mechanism underlying evolu- 
tion in ecological systems, but we are often unable to 
ans"wer questions concerning stability and the nature of 
the dynamical evolution (intermittent versus gradual) . It 
is also difficult to measure the degree of interrelatedness 
of an ecology (e.g. Bj0rnstad et al. , 2001): "who is in- 
teracting "with "whom and ho"w strongly, and it is difhcult 
to determine the species abundance in detail. What is 
especially difficult is to monitor the temporal variation 
in the aforementioned quantities over evolutionary time. 

Can general principles be identified for the overall dy- 
namical behaviour of evolution? Even if the character- 
istics of each individual species have to be considered in 
their proper specific contexts, perhaps general la"ws do 
operate at an overall level. Obviously the ans"wers to 
these questions are empirical, but indicators may be ob- 
tainable from deliberately simplified theoretical models. 
It is obviously important to consider carefully the type of 
simplification assumed. Simplified models often operate 
directly at species level and typically consider only a fe"w 
coupled species (see e.g. Doncaster, 2000) but in order 
to capture the consequences of the complexity character- 
istic of ecology "we believe it is important to treat species 
as emergent structures and to allo"w for the multitude 
of interactions each individual (and therefore species) is 
subject to. We find that the very nature of the dynamics 
of the ecology is strongly influenced by the complexity of 
the system. 

Temporal as "well as spatial variations in the physical 
environment are kno"wn to play an important role in evo- 
lution. It is also often assumed that co-evolution "with in- 
teraction bet"ween co-existing individuals, or species, may 
influence the evolutionary dynamics in a significant "way 



(Kauffman, 1995; Bak & Sneppen, 1993). The relative 
importance of selective pressure of purely physical origins 
and co-evolutionary effects is not clear, and it seems dif- 
ficult to resolve the issue solely by use of selected specific 
case studies. Moreover, seen from an ecological point of 
vie"w, the biotic and the physical environment are cou- 
pled. 

In the present paper "we present an individual based 
mathematical model of an evolving ecology. The model 
is kept sufficiently simple to allo"w us to simulate evolu- 
tionary time scales. 

We attempt here to gain some insight into the possi- 
ble effects of co-evolution through the study of a model 
in "which variations in the physical environment are alto- 
gether neglected. Our model is not meant to be a realistic 
representation of biological evolution, but rather a theo- 
retical approximation in which co-evolution is made to be 
the prominent driving force. We then demonstrate within 
this model that speciation does occur and we study in 
some detail the dynamical features of the evolution of 
the model as well as the nature of the ecology created by 
the co-evolutionary dynamics. 

We are interested in the qualitative behaviour of a sys- 
tem in which the mutual interaction between co-existing 
individuals of different genetic composition determines 
the possibility of the individual to thrive. The model 
emphasises the web of interactions between individuals 
of different genomic composition, to stress this aspect 
we will talk about the Tangled Nature model - or the 
TaNa model for short. We represent biotic factors in 
terms of the co-evolutionary effects on the fitness of in- 
dividuals. The model is a simplification. No distinction 
is made between genotype and phenotype and the de- 
tails of the reproductive mechanism are kept to a min- 
imum. This simplification allows us to represent evolu- 
tion in terms of the dynamics of the distribution of the 
population in genome space. We demonstrate that at 
a qualitative level the complex dynamics of the model 
resembles known aspects of long term biological evolu- 
tion such as speciation and intermittent behaviour. We 
are also able to study the competition between asexual 
and sexual reproduction. We find that asexual repro- 
duction is superior during periods of rapid relocation of 
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the configuration in genome space, whereas sexual repro- 
duction is most advantageous during the coherent more 
stable epochs. It is natural to relate these stable epochs 
to periods of Evolutionary Stable Strategies (ESS), as in- 
troduced by Maynard Smith (1982). The stable periods 
of our model are, however, not perfectly stable as fluc- 
tuations caused by mutations can trigger a switch from 
one stable period to another. We therefore suggest call- 
ing these periods "quasi-Evolutionary Stable Strategies" 
or q-ESS. The overall effect of the evolutionary dynamics 
of the present model is to increase the average duration 
of the q-ESS. 

II. DEFINITION OF THE MODEL 
Interaction 

We now define the TaNa model in detail. We repre- 
sent individuals in the same way as in models consid- 
ered by, e.g., Kauffman (1995), Higgs and Derrida (1992), 
Gavrilets (1999), Eigen et al. (1988), and Wagner et 
al. (1998). An individual is represented by a vec- 
tor S" = (S'f,Sf,...,S'£) in genotype space. Here Sf 
may take the values ±1. These may be interpreted as 
genes with two alleles, or a string of either pyramidines 
or purines. Individuals are labelled by Greek letters 
Q!,/9, ... — 1,2, N{t). When we refer, without refer- 
ence to a specific individual, to one of the 2^ positions 
in genome space, we use roman superscripts S°, S**, ... 
with a, 6, ... = 1,2, ...,2^. Note, many different individ- 
uals S",S'^,. . ., may reside on the same position, say 
in genome space. Geometrically the vector S" repre- 
sents one of the corners of the L dimensional hyper-cube 
5 = {-1,1}^ (see Fig. 1). 




FIG. 1. A three dimensional genome space. For L = 3 the 
sequence of genes uniquely defines a vertex of a cube. The 
number of edges (dotted lines) that must be traversed be- 
tween 2 vertices defines their Hamming distance. Interactions 
between vertexes are shown as solid curves with thickness in- 
dicating the strength and circles placed at the vertexes have 
radii proportional to the occupation (number of individuals 
present) with the genome in question. Note that interactions 
are defined even for unoccupied vertices. 

The ability of an individual a to reproduce is controlled 
byi?(S",t): 

^(S",t) = ^ E E ^.(S",S'^)5f 5f - f,Nit), (1) 

where N{t) is the total number of individuals at time t. 
The sum over individuals /3 in Eq. (|^) is more conve- 
niently expressed as a sum over the locations S in the 
genome space S, using the occupancy n{S, t) of the loca- 
tions we obtain: 

1 ^ 

i?(S",t) = ^ ^ MS'^,s)srs, - f^N{t). 

(2) 

Two positions S° and S** in genome space are coupled 
with the fixed strength J"*" = J(S°,S''). This coupling 
is non-zero with probability 0, in which case we assume 
jafc _^ j6a -j-Q |-jg random and uniform on the interval 
[— c, c], where c is a constant. The structure of the cou- 
pling in genome space is sketched in Fig. 1. 

Some comments about the interaction matrix J(S°, S**) 
are appropriate. In our simplistic approach a given 
genome is imagined to lead uniquely to a certain set of at- 
tributes (phenotype) of the individuals/organisms. The 
locations S° and represent blueprints for organisms 
that exist in potentia. The locations may very likely be 
unoccupied but, if we were to construct individuals ac- 
cording to the sequences S" and S** the two individuals 
would have some specific features. Anecdotally we can 
imagine that S° corresponds to rabbits and S'' represents 
foxes. The number J(S'',S'') now represents the poten- 
tial influence of an individual constructed according to 
the genome sequence S** on an individual constructed ac- 
cording to the genome sequence S". In our toy example 
J(S°, S*) represents the fact that the foxes will tend to 
eat the rabbits and thereby decrease the rabbits ability to 
survive and J(S'',S°) represents the fact that the avail- 
ability of rabbits as a food source will help to sustain the 
foxes. Other examples could be parasitic or collabora- 
tive relationships. In order to emphasis co-evolutionary 
aspects we have excluded "self-interaction" among indi- 
viduals located at the same positions S in genome space, 
i.e., we use J(S, S) = 0. 

In reality the mutual influence between two individuals 
of a certain genotype (phenotype) is, of course, not a ran- 
dom quantity. The interaction maybe be collaborative. 
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competitive or neutral. It is this aspect we represent 
by ascribing a set of fixed randomly assigned coupling 
strengths between the positions in genome space. 

We stress that the segregation (or speciation) to be 
discussed below is an effect of different couplings be- 
tween different positions S° and S**. When we assume 
J,(S'^,S'') = Jo independent of and S^ we find the 
population not to be concentrated around a subset of 
all positions in genome space, instead the population is 
smeared out through the genome space in a diffuse man- 
ner. 

The conditions of the physical environment are simplis- 
tically described by the term nN{t) in Eq. (1), where 
determines the average sustainable total population size. 
An increase in /i corresponds to more harsh physical con- 
ditions. This is a simplification, though one should re- 
member that the physical environment encountered by 
an organism is to some extent produced by the presence 
of other living organisms. Consider, for example, the 
environment experienced by the bacterial flora in the in- 
testines. Here one type of bacteria very much live in an 
environment strongly influenced by the presence of other 
types of bacteria. In this sense some fluctuations in the 
environment may be thought of as included in the matrix 
J(S°,S''). 
Reproduction 

Asexual reproduction consists of one individual being re- 
placed by two copies, this event occurs for individuals S" 
with a probability per time unit proportional to 



Po//(S",t) 



exp[f/"(S",i)] 
1 + exp[i?(S",i)] 



e [0, f ] 



(3) 



In the case of sexual reproduction an individual S" 
is picked at random and paired with another ran- 
domly chosen individual with Hamming distance d 

2 



IS"" - Sh < d 



2 ^1=1 '^i \ - "^max (allowing at most d,nax pairs 

of genes to differ) . The pair produces an offspring 7 with 
a probability ^po//(S", t)po//(S'^, t), with S^ chosen at 
random from one of the two parent genes, either S" or 

Mutation 

Genes mutate with probability Pmut, represented by a 
change of sign —S], during the reproduction pro- 

cess. Choosing genes at random from the parents may 
be thought of as a process similar to recombination for 

dfnax ^ 2. 

Annihilation 

For simplicity an individual is removed from the system 
with a constant probability pkui per time step. This pro- 
cedure is implemented both for asexual and sexually re- 
producing individuals. 
Time Step 

A time step consist of one annihilation attempt followed 
by one reproduction attempt. One generation consists of 
N{t)/pkiu time steps, the average time taken to kill all 
currently living individuals. 



Stability 

At an average level of description, and neglecting muta- 
tions, the above dynamics is described by the following 
set of equations (one equation for each position in the 
genotype space): 



9n(S,t) 
dt 



[po//(S,t) -pkiii]n{S,t) 



(4) 



controlling the temporal evolution of the occupancy 
n{S,t) of the positions S in genotype space S. Station- 
ary solutions (i.e., those for which dn/dt = 0) demand 
either n{S,t) = or Poff{S,t) = Pkui- During the q-ESS 
the system manage to find a configuration in genotype 
space for which all occupied positions satisfy the balance 
between production of offspring and decease. The fit- 
ness poff{S"-,t) of individuals at a position S° depends 
on the occupancy n{S'',t) of all the sites S^ with which 
site S° is connected through couplings J"''. Accordingly, 
a small perturbation in the occupancy at one position is 
able to disturb the balance between poff{S,t) and pkui 
on connected sites. In this way an imbalance at one site 
can spread as a chain reaction through the system, pos- 
sibly affecting a global reconfiguration of the genotypical 
composition of the population. 



III. DYNAMICAL BEHAVIOUR 

We consider three different types of populations. (1) 
a purely asexual population, (2)purely sexual population 
and (3) a mixed population in which mutations can trans- 
form an asexual individual into a sexually reproducing 
individual and vice versa. 



A. Asexual reproduction 

In this subsection we discuss the model when all repro- 
duction is assumed to be asexual. 
Initiation 

Let us consider the initiation of the model. First we 
place the entire population A^(0) at a randomly chosen 
location S* in genome space. The ff-function in Eq. (|^) 
wiU be given by H{S*,0) = -yuiV(O) since n(S) = for 
S ^ S* and J(S*, S*) = 0. If no mutations can occur the 
population will remain confined at the location S* and 
the size of the population n{S*,t) will according to Eq. 
(0) approach the value 



N* = ih/^"^^'" 



PkiU 



Mutations do occur, however, and the population will 
migrate away from the original location S* into the sur- 
rounding region of genome space. In Fig. 2, we show a 
cladogram indicating the evolution of the first 110 gen- 
erations. During this initial period the newly invaded 
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positions are only occupied for a few generations. After 
this period of rapid changes a relatively stable configura- 
tion is achieved, and the occupied positions to the right 
in Fig. 2 indicate that the system has entered its first 
q-ESS. 



Root- 




New q-ESS is 
established. 
d„-~9-10 
from root position 



500 Identical 
Individuals 



Time 



FIG. 2. The initial diversification from a single position in 
genome space. The system is initialised at time t = with 500 
identical individuals and allowed to develop autonomously. 
The system to mutates away from the initial location, which 
becomes extinct relatively quickly. After 34 branchings the 
system finds a stable configuration and enters the first q-ESS 
(see Figs. 3 and 4). 

We have also studied simulations started out from an 
initial population spread out over many randomly chosen 
positions in genome space. Most of these initially occu- 
pied positions rapidly become extinct. In this way, the 
diversity in genome space passes through a "bottleneck" 
before the population starts to migrate out into genome 
space from one or a few positions which were able to pass 
through the bottleneck. Prom then on, the evolution of 
the ecology behaves in the same way as when started out 
from one single position in genome space. 
Long time behaviour 

Now we turn to a discussion of the nature of the long 
time dynamics of the model. The model consists of a 
variable number of co-evolving individuals all subject to 
the same physical environment. An individual's ability to 
thrive depends on its own genetic composition as well as 
the genetic composition of the other individuals present. 
The dynamical evolution, driven by mutations, will have 
to strike a balance between the multiplication of the in- 
dividuals and the total carrying capacity of the environ- 
ment. Different types of genotypical compositions of the 



population can achieve this balance. 

One possibility consists of very numerous populations 
distributed on a relatively small number of isolated re- 
gions in genotype space corresponding to a small num- 
ber of species (compare to the total number of genotypes 
for a given genome length). These configurations can be 
stable for very many generations and allow the species 
to co-exist quietly during coherent periods of little vari- 
ation in the total size or composition of the population, 
see Fig. 3. 




100000 200000 300000 400000 500000 

No. Generations 
FIG. 3. Change of total population with time for a system 

with L = 20, H = 0.1375, Pk^ll = 0.2, Pmutate = 0.01, c = 100 

and = 0.25. Regions of high population and low relative 
fluctuations (q-ESS) are clearly distinct from regions of low 
population and high relative fluctuations (transition periods) . 

In Fig. 4, we demonstrate that the occupancy of the 
positions in genome space fluctuates only very little dur- 
ing these stable periods. We call these epochs q-ESS, or 
quasi Evolutionary Stable Strategies (Maynard Smith, 
1982). The q-ESS exhibit a degree of stability against 
mutation induced changes, but fluctuations in the fre- 
quency distribution in genome space can abruptly desta- 
bilize such a configuration. 
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100000 200000 300000 400000 500000 
No. Generations 
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FIG. 4. Occupation of genome space versus time for the 
same simulation as in Fig. 3. We arrange the positions in 
genome space in a convenient arbitrary way along the y-axis 
and place a dot for each occupied location at a given time. 
Periods of stability (q-ESS) interrupted by periods of hectic 
rearrangement are clearly visible. 

Wc show, however, in Fig. 5 that the distribution of 
lifetimes of the q-ESS, measured in numbers of genera- 
tions, is very broad. 
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FIG. 5. Log-Log plot of the distributions of lengths of 
q-ESS (solid line), transition periods (dotted lino) and life- 
times of occupied locations in genome space (dot-daslicd line) . 
We observe power-law-like behaviour in both the q-ESS and 
the lifetimes of genome space locations, but the transition pe- 
riods exhibit an abrupt cut-off at much shorter times. The 
lifetimes curve extends further than that of the q-ESS, indi- 
cating that locations may remain occupied from one q-ESS to 
another, surviving the transition. 

Transitions 

The q-ESS periods arc separated by periods of hectic 
rearrangements of the geiiotypical composition of the en- 
tire population. During these periods of rapid change, 
the total number of individuals is small and populations 
located at specific positions in genome space undergo se- 
quences of bifurcations as seen in Fig. 6, where we follow 
the evolution across the hectic transition period from one 
q-ESS to the next. 



Root 





1660 1690 
Time 

FIG. 6. The descendants of a fit location across a transi- 
tion period. The Lower portion of this figure is a slice from an 
occupation plot, similar to Fig. 4. Wc track the descendants 
of a single, fit location across the transiton and into the next 
q-ESS. After 12 Branchings the descendants have found a new 
fit configuration that is stable enough to form part of the new 
q-ESS. The original location does not survive the transition 
and lineages of the other fit locations from the original state 
die out very rapidly. 

The figure is a cladogram tracing out all the descen- 
dants originating from one root. One notices that most 
of the new positions spun off from the root die before 
the next q-ESS is reached. While new branches are cre- 
ated old ones die. The periods of rapid rearrangement 
in genome space are transition periods during which the 
system searches for a new stable configuration. 
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FIG. 7. Flower diagrams of the configuration of a sys- 
tem during a transition from one coherent state to another. 

Flower diagrams visualise the interactions and genome space 
proximities in a system at a given time. Very fit locations 
have a complete set of nearest mutational neighbours (these 
make up a "flower"). Each occupied location is represented 
by a circle of radius proportional to its occupation, (or a num- 
ber, for very large locations) positive interaction strengths are 
solid lines, and negative interaction strengths dotted lines. If 
a flower is out competed by a new mutant, the q-ESS is dis- 
turbed and the system executes an adaptive walk, searching 
for a new q-ESS. Each diagram is separated by approximately 
5 to 10 generations. See the main text for discussion. The 
Hamming distance between two highly occupied positions ap- 
pears midway between them. 

The nature of the transition from one q-ESS to the 
next is indicated in Fig. 7. This set of diagrams repre- 
sent in a quantitative way the positions with the largest 
occupation together with the couplings in genome space. 

Gene sequences in the diagrams are as follows: 



8^ = 


(+1, +!,+!, 


-1, 




82 = 


(+1,+1,-1, 


-1, 




83 = 




-1, 


-1,+1, 


8^ = 




-1, 


-1,+1, 


8^ = 




-1, 




8^ = 


(+!,+!, +1, 


-1, 


-1,+1, 


8^ = 


(+1,+1,-1, 


-1, 


-1,+1, 


8^ = 




-1, 


-1,+1, 


S« = 


(-1.+1.+1. 


-1. 


-1,+1, 


S"' = 


(-1,-1,-1, 


-1, 


-1,+1, 


8^1 = 




-1, 




8^2 = 


(-!,+!, -1,+1, 





During a transition between one q-ESS and the next, 

the systems behaviour becomes very hectic. Starting at 
(a) we see that a new mutant has invaded the previous 
coherent configuration (originally similar to (b)) with a 
negative interaction with most of the existing flowers but 
a strong enough positive interaction with one of them for 
it to survive. This causes the coherent state to be desta- 
bilised. In (b) we see that although the new mutant does 
not survive for long, it has drastically reduced the popu- 
lation at 8^, which in turn has a harmful effect on 8^ and 
8^. In (c) we observe that two further new mutants have 
been able to invade, this is due to the reduced fitness of 
the original sites from the effect of the first invader. The 
new mutants are transient, they represent steps on an 
adaptive walk. The system is now in a situation where 
it is partly executing such a walk, and partly still in the 
previous coherent state. This continues into (d), where 
we can see that 8^ and 8"^ are still holding on, and their 
complete first circles evince they are still reproducing. 
By (e), however, things have changed again. The adap- 
tive walkers are now out-competing the originals, 8^ has 



become extinct and 8"^ has a very low population. We 
also observe the formation of a double flower (8^ and 
8*) which consists of two fit centres in adjacent locations 
in genome space. The adaptive walk continues for some 
time until a new q-ESS is found at (f). 

We have studied the distribution of non-zero couplings 
J(8'',8'') between a given occupied position 8" and an- 
other occupied position 8*'. During the q-ESS this distri- 
bution is narrow and its average is smaller than during 
the transition periods, where the distribution broadens. 
Epoch Distributions 

It is interesting to take a further look at Fig. 5. One 
notices that the distribution of lifetimes of occupied po- 
sitions reaches as far out as the distribution of q-ESS 
durations. In fact we observe in the simulations that po- 
sitions sometimes are able to remain occupied across the 
transition from one q-ESS to the following, correspond- 
ing to a species that survives a mass extinction. Fig. 5 
also shows clearly that the periods of hectic reconfigu- 
ration typically last for a significantly smaller number 



of generations than do the q-ESS periods. Finally, it is 
very interesting to note that both the lifetimes of individ- 
ual positions and the distribution of q-ESS epoch lengths 
are power-law-like with exponents around -2.3 and -1.8 
respectively. We mention that the distribution of q-ESS 
durations can be compared to the distribution of lifetimes 
of genera obtained from the fossil record. The latter has 
a shape similar to the distribution of q-ESS durations 
shown in Fig. 5. Power-law fits to the fossil record data 
leads to an exponent around 2. For a recent analysis of 
data from the fossil record see Newman & Sibani (1999). 
Adaptation Level Increases 

We now turn to a discussion of the overall long time ef- 
fect of the dynamics of the TaNa model. How does the 
genomic composition of early configurations differ from 
those generated after hundreds of thousands of genera- 
tions? 
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FIG. 8. Running Averages of the duration of q-ESS and 
transition periods. At the end of a particular state, we eval- 
uate the average length of states up until that point. We see 
that the average length of the transition periods settles down 
and fluctuates slightly around constant, whereas the average 
q-ESS length continues to increase. 

In Fig. 8 we show the running average of the durations 

of the q-ESS as well as the transition epochs. One notices 
that there is no significant trend in the duration of the 
hectic periods of rearrangement separating the consecu- 
tive q-ESS. The average duration of the q-ESS periods, 
however, slowly increases with time. This means that the 
entire ecology gradually becomes more stable. Or we may 
say that the ecology (represented by the distribution of 
the population through genome space) becomes increas- 
ingly better adapted; not adapted to some fixed external 
environment, but adapted in the sense that the ecology 
as a whole achieves collectively increasingly stable con- 
figurations among the total set of all possible ways of dis- 
tributing a population through genome space. Does this 
mean that eventually some maximally "fit" or adapted 
configuration is reached? Our simulations indicate, as 
expected, that the time to reach a stationary state in- 
creases exponentially with increasing genome length L. 
We will accordingly expect that for biologically relevant 
systems an ecology would never have the time to reach 
a final stationary state. Moreover, even if the system 
becomes stationary in the sense that the average dura- 
tion of the q-ESS becomes time independent, switching 
between different equally well adapted configurations is 
likely to continue forever. From the statistical mechan- 
ics of disordered systems we do not expect the optimally 
adapted configuration to be unique. Hence transitions 
between equally maximally adapted configurations may 
continue even in the mathematical limit of infinitely long 
time. 

The increase of the average duration of the q-ESS can 
be viewed as an optimising process. This is in accor- 
dance with the suggestion (Mayr, 1988) that the effect 
of biological evolution is to optimise some quantity. The 
identification of the quantity being optimised is still de- 
bated (Fogel & Beyer, 2000). Unfortunately we cannot 
identify a specific mathematical function of the distribu- 
tion n{S,t) in genome space which is optimised as an 



effect of the dynamics. However, it is very interesting to 
relate the average duration of the q-ESS to the extinc- 
tion rate. Due to insufficient statistics we cannot, un- 
fortunately, make a quantitative comparison. We note, 
qualitatively, that an increasing average duration of the 
q-ESS corresponds to a decreasing extinction rate. This 
is consistent with Raup & Sepkoski's (1982) analysis of 
the fossil record, which suggests that the extinction rate 
might have declined through the Phanerozoic. 



B. Sexual Reproduction 

We now briefly discuss a model in which all individuals 
are assumed to reproduce sexually. More detail will be 
the presented in a future communication. 
Long Time Behaviour 

In Fig. 9a and Fig. 9b, we show the temporal behaviour 
of the total number of individuals together with the oc- 
cupancy in genome space. 
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FIG. 9. The total number of individuals (below) and the 
occupation plot of the species (above). The horizontal axis is 
the generational time. The different plateaus in the popula- 
tion size correspond to rearrangements of the population in 
genome space. Parameters arc L — 20, jj, = 0.07, pkUl = 0.2, 
Pmutate = 0.01, c = 100 and e = 0.25. 

We have assumed that the maximum number of genes 
which parents can differ, dmax = 2 and in Fig. 9b, we 
plot only the species occupancy, that is we have coarse 
grained genome space with a resolution of dmax- This is 
done in the following way: In each time step we identify 



8 



the position with the largest population, we lump this po- 
sition together with all positions within a distance dmax ■ 
Next we find the position with the second largest popu- 
lation, and lump this location together with all positions 
within distance dmax- We continue this until all occupied 
positions have been considered. The locations in genome 
space are labelled in a convenient, but arbitrary, way. For 
each time step we place a dot along the y-axis for each 
occupied (coarse grained) positions in genome space. Fi- 
nally, along the x-axis, we convert time steps into time 
measured in generations. We observe that, similar to the 
asexual case, the model evolves through a set of q-ESS 
phases separated by short transition periods. We also 
emphasise that well established species can be identified 
as the well separated locations in genome space where 
the population is concentrated. 
Lifetime Statistics 

In Fig. 10, we show the distribution of lifetimes of oc- 
cupations of individual multiple occupied positions in 
genome space. A slow power-law-like decay is observed. 
Note the similarity with the distribution found in the 
asexual case and with the distributions reported from 
the fossil record , see e.g. Newman & Sibani (1999). 
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FIG. 10. Distribution of positions's lifetime (log-log plot) 
for the sexual case. Same simulation as in Fig. 9. 

C. Competition Between Sexual and Asexual 
Populations 



spring. This is compensated by assuming the asexual 
mutation rate for the reproductive mode gene to be half 
that for the sexual. In this way we eliminate any net drift 
induced by the rates of mutation from one reproductive 
mode to the other. 

In Fig. 11, we plot the total population of the sys- 
tem, along with those of the two subpopulations (i.e. the 
numbers of sexual and asexual individuals present at a 
given time) we see that in the coherent phases the sys- 
tem is predominantly asexual, despite the large fluctua- 
tions in population, whereas during coherent phases we 
see the opposite: the population becomes predominantly 
sexual. The reason for this effect is not clear, but we 
believe that the tight structures (clusters of neighbour or 
nearest neighbour positions all of significant occupancy) 
evolved in genome space by a sexual population may be 
more suited to a constant environment than the more 
scattered ones observed for asexual populations. 

This hypothesis is supported by the observation that a 
sexual population with dmax — L shows q-ESS with asex- 
ual dominance (see Fig. lib). A sexual population of this 
sort would not necessarily form such tight stable struc- 
tures since the pressure to have nearby mate-genomes of 
high fitness would not be present. 

The observed asymmetry between the asexual and the 
sexually reproducing subpopulations for small values of 
dmax may also be related time scales. The sexually re- 
producing individuals need other individuals in the im- 
mediate vicinity of genome space which takes time to es- 
tablish. However, the mixing of "genes" involved in the 
recombination during sexual evolution allows the sexual 
population to scan through larger portions of genome 
space, and therefore, if given time, this population is 
likely to find the better adapted configurations. 
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To observe more directly the differences in behaviour 
between sexual and asexual populations we have con- 
structed a mixed reproductive mode model. Here, indi- 
viduals are given an extra gene which does not explicitly 
enter the Hamiltonian, but instead dictates an individ- 
ual's reproductive mode. Mutations to this gene occur 
during reproduction in the normal way. An asexual par- 
ent may potentially produce two sexual offspring whereas 
a sexual parent may produce at most one asexual off- 
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FIG. 11. Population Sizes versus Time 

in Mixed-Reproductive-Mode Systems. Again we see a di- 
vision in behaviours between q-ESS and transition periods, 
(a) During q-ESS the sexual population dominates, whereas 
the reverse is observed during transitions. Evidence to sup- 
port the hypothesis that this is due to tightly packed flower 
structures evolved by a sexual population is supported by (b) 
in which the selection pressure to form these structures is re- 
moved by allowing all sexual individuals to reproduce. The 
arrow marks a q-ESS in which the asexual population domi- 
nate. Others q-ESS where the asexual population dominates 
are observed on shorter timescales than are visible on this 
plot. 

IV. DISCUSSION AND CONCLUSION 

We have discussed a simple, very general mathematical 
metaphor (Gravrilets 1999) by which we can study the 
long time (of order 10^ or 10^ generations) behaviour of 
an ecology. 

Both asexual and sexually reproducing populations 
evolve through a set of relatively stable configurations, 
the q-ESS, separated by short transition periods of hectic 
reorganisation of the genomic composition of the ecology. 
The population segregates in genome space into well sep- 
arated clusters of highly occupied positions. Speciation 
events occur when a position or a tight cluster of posi- 
tions undergoes successive bifurcations in genome space. 
This type of behaviour is observed for a broad range of 
control parameters. 

The co-evolutionary dynamics produce a highly tan- 
gled interdependent population of species. The evolution 
gradually increases the robustness of the entire ecology 
against fluctuations in the genomic and physical envi- 
ronment. In agreement with analysis of the fossil record 
we find that the average duration of the q-ESS increases 
slowly with time. 
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